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SPECTRAL ANALYSIS OF THE TRUNCATED HILBERT 
TRANSFORM WITH OVERLAP 



REEMA AL-AIFARli AND ALEXANDER KATSEVICH^ 



Abstract. We study a restriction of the Hilbert transform as an operator Hj- 
from L^{a2, 04) to L'^{ai, 03) for real numbers ai < a2 < < a^. The opera- 
tor Hj' arises in tomographic reconstruction from Umited data, more precisely 
in the method of differentiated back-projection (DBF). There, the reconstruc- 
tion requires recovering a family of one-dimensional functions / supported on 
compact intervals [02, a4] from its Hilbert transform measured on intervals 
[01,03] that might only overlap, but not cover [02,04]. We show that the in- 
version of Ht is ill-posed, which is why we investigate the spectral properties 
of Ht- 

We relate the operator Ht to a self-adjoint two-interval Sturm- Liouville 
problem, for which we prove that the spectrum is discrete. The Sturm-Liouville 
operator is found to commute with Ht, which then implies that the spectrum 
of H^Ht is discrete. Furthermore, we express the singular value decompo- 
sition of Ht in terms of the solutions to the Sturm-Liouville problem. The 
singular values of Ht accumulate at both and 1, implying that Ht is not 
a compact operator. We conclude by illustrating the properties obtained for 
Ht numerically. 



1. Introduction 

In tomographic imaging, which is widely used for medical applications, a 2D 
or 3D object is illuminated by a penetrating beam (usually X-rays) from multiple 
directions, and the projections of the object are recorded by a detector. Then one 
seeks to reconstruct the full 2D or 3D structure from this collection of projections. 
When the beams are sufficiently wide to fully embrace the object and when the 
beams from a sufficiently dense set of directions around the object can be used, 
this problem and its solution are well understood [T^]- When the data are more 
limited, e.g. when only a reduced range of directions can be used or only a part 
of the object can be illuminated, the image reconstruction problem becomes much 
more challenging. 

Reconstruction from limited data requires the identification of specific subsets of 
line integrals that allow for an exact and stable reconstruction. One class of such 
configurations that have already been identified, relies on the reduction of the 2D 
and 3D reconstruction problem to a family of ID problems. The Radon transform 
can be related to the ID Hilbert transform along certain lines by differentiation 
and back- projection of the Radon transform data [differentiated hack-projection or 
DBF). Inversion of the Hilbert transform along a family of lines covering a sub- 
region of the object {region of interest or ROI) then allows for the reconstruction 
within the ROI. 
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This method goes baek to a resuh by Gelfand and Graev Its apphcation 
to tomography was formulated by Finch [4] and was later made explicit for 2D in 
[H [H US] and for 3D in [H [Til [H [M] . To reconstruct from data obtained by 
the DBF method, it is necessary to solve a family of ID problems which consist of 
inverting the Hilbert transform data on a finite segment of the line. If the Hilbert 
transform H f of a ID function / was given on all of M, then the inversion would be 
trivial, since H^^ — —H. In case / is compactly supported, it can be reconstructed 
even if H f is not known on all of M. Due to an explicit reconstruction formula by 
Tricomi pTj , / can be found from measuring H f only on an interval that covers the 
support of /. However, a limited field of view might result in configurations in which 
the Hilbert transform is known only on a segment that does not completely cover 
the object support. One example of such a configuration is known as the interior 
problem [T|[9l lllll20) . Given real numbers ai < 02 < as < 04, the interior problem 
corresponds to the case in which the Hilbert transform of a function supported on 
[01,04] is measured on the smaller interval [02,03]. 

In this paper, we study a different configuration, namely supp / = [02,04] and 
the Hilbert transform is measured on [01,03]. We will refer to this configuration 
as the truncated problem with overlap: the operator Ht we consider is given by 
'P[ai,a3]HV[a2.ai]i whcrc H is the usual Hilbert transform acting on and 
Vn stands for the projection operator {Vnf){x) = f{x) if x £ ft, {Vnf){x) — 
otherwise. (The interior problem corresponds to f2i C ^2', one can also consider 
operators of the type VniHVn2 where fii and ^2 are separated by a gap, as in 
[7j; we would call these truncated Hilbert transform with a gap). Figure [T] shows 
the different setups. Examples of configurations in which the truncated Hilbert 
transform with overlap and the interior problem occur are given in Figures [2] and 
[21 The truncated problem with overlap arises for example in the "missing arm" 
problem. This is the case where the field of view is large enough to measure the 
thorax but not the arms. 

suppf suppf 



b a, 

Hf measured Hf measured 

b 33 
(a) Finite Hilbert transform. (b) Interior problem. 



suppf suppf 

32 34 3^ 32 

Hf measured Hf measured 



3^ 33 33 34 

(e) Truncated Hilbert transform with overlap, (d) Truncated Hilbert transform with a gap. 



Figure 1. Different setups for VniHVn2- The upper interval 
shows the support 0,2 of the function / to be reconstructed. The 
lower interval is the interval ili where measurements of the Hilbert 
transform Hf are taken. This paper investigates case (c). 
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Figure 2. Two examples of the truncated problem with overlap, 
Fig. 2(a) shows the missing arm problem. In both cases, the 



field of view (FOV) does not cover the object support. On the line 
intersecting the object, measurements can only be taken within the 
FOV, i.e. from oi to 03. The Hilbert transform is not measured 
on [as, 04]. Consequently, a reconstruction can only be aimed at in 
the grey-shaded intersection of the FOV with the object support, 
called the region of interest (ROI) . 




Figure 3. The interior problem. Here, the FOV also does not 
cover the object support. The line intersecting the object is such 
that the Hilbert transform is only measured in a subinterval [02, 03] 
of the intersection [oi, 04] of the line with the object support. The 
ROI is the grey-shaded intersection of the FOV with the object 
support. In this case stable reconstruction of the shaded ROI is 
impossible unless additional information is available. 
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Fix any four real numbers oi < 02 < 03 < 04. We define tlie truncated Hilbert 
transform with overlap as the operator 

Definition 1. 

{HTf){x) := -p.v. / ^^dy, a; €(01,03) (1-1) 

7^ Ja2 

where p.v. stands for the principal value. In short, 

Ht ■■= 'P[ai^a3]HV[a2,a4], 

where H is the ordinary Hilbert transform on (R) . 

As we will prove in what follows, the inversion of Ht is an ill-posed problem 
in the sense of Hadamard [3]. In order to find suitable regularization methods 
for its inversion, it is crucial to study the nature of the underlying ill-posedness, 
and therefore the spectrum a{H^HT). An important question that arises here 
is whether the spectrum is purely discrete. This question has been answered for 
similar operators before, but with two very different answers. In |10| . it was shown 
that the finite Hilbert transform defined as Hp = 'P[a,b]HV[a,b] has a continuous 
spectrum a{HF) = [—i, i]- On the other hand, in [8], we find the result that for the 
interior problem Hj = 'P[a2,a3]HV[ai,a4], the spectrum a{HjHi) is purely discrete. 

The main result of this paper is that H^Ht has only discrete spectrum. In 
addition, we obtain that and 1 are accumulation points of the spectrum. Fur- 
thermore, we find that the singular value decomposition (SVD) of the operator Ht 
can be related to the solutions of a Sturm-Liouville (S-L) problem. For the actual 
reconstruction, one would aim at finding / in (jl.ip only within a region of interest 
(ROI), i.e. on [02,03]. A stability estimate as well as a uniqueness result for this 
setup were obtained by Defrise et al in [21 . A possible method for ROI reconstruc- 
tion is the truncated SVD. Thus, it is of interest to study the SVD of Ht also for 
the development of reconstruction algorithms. 

In [7] and [5], singular value decompositions are obtained for the truncated 
Hilbert transform with a gap 'P[a3.a4]-ff^[ai,a2] ^-^id fo^' ^i- This is done by re- 
lating the Hilbert transforms to differential operators that have discrete spectra. 
We follow this procedure, but obtain a differential operator that is different in na- 
ture. In [7| and [8] the discreteness of the spectra follows from standard results of 
singular S-L theory (see e.g. [H]). In the case of truncated Hilbert transform (jl.ip 
we have to investigate the discreteness of the spectrum of the related differential 
operator explicitly. 

The idea is to find a differential operator for which the eigenfunctions are the 
singular functions of Ht on (02,04). We define the differential operator similarly 
to the one in , |H] , but then the question is which boundary conditions to choose 
in order to relate the differential operator to Ht- To get an intuition about the 
singular functions of Ht we use the well-known fact that the Hilbert transform 
of the characteristic function X[a,b] has logarithmic singularities at the points of 
discontinuity: 

/■ 1 X[a,6](j/) , 1, 

p.v. / ■ — ■ ay = — In 

J TT y — X TT 

Let {(7nT fn, gn} denotc the singular system of Ht that we want to find. The 
problem can be formulated as finding a complete orthonormal system {/rajnsN in 
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£^(02,04) and an orthonormal system {gnjneN in £^(01,03) such that there exist 
real numbers Un f^or which 

Hrfn — ^n9n-, 
Hj'g-n — CTnfn- 

At the moment, the g^s only have to be complete in Ran(i7T), but as we will see 
in Section[5l Ran(i/T) is dense in £^(01,03). 

As will be shown in Sectional the functions /„ and .g„ can only be bounded or of 
logarithmic singularity at the points . Suppose /„ has a logarithmic singularity at 
aj. Since Ht integrates over [02, 04], the function Ht/v. would be neither bounded 
nor of logarithmic singularity at 02- Hence, this would violate the property of gn 
at a2- Therefore, /„ has to be bounded at 02- If fn does not vanish at aj, this 
leads to logarithmic singularities of HtJu and g„ at 02. Using the same argument 
we conclude that g„ is bounded at and /„ has a logarithmic singularity at 03. 

On the other hand, since g„ is bounded at 03, Hxfn is also bounded there. This 
requires that close to 03, /„ ~ fn.i + /n.2ln|a:; — a^l for functions fn.i continuous 
at 03. A similar argument holds for g„ at 02- Close to that point, gn = gn.i + 
gn,2 In |a; — 02] for functions gn,i continuous at 02. 

Clearly, HtJu is bounded at ai and Hxgn is bounded at 04. Therefore, /„ has 
to be bounded at 04 and gn must be bounded at ai . 

Thus, if we want to show the commutation of Ht with a differential operator 
that acts on fn{x), x e (02, 04), we need to impose boundary conditions at a2 and 
04 that require boundedness and some transmission conditions at 03 that make the 
bounded term and the term in front of the logarithm in /„ continuous at 03. 

bdd bdd 




(a) Sketch of fnS. (b) Sketch of gn's. 



Figure 4. Intuition about the singular functions of Ht- 



Having found these properties of the singular functions of Ht (in case the SVD 
for Ht exists), in Section [5] we introduce a differential operator and find a self- 
adjoint extension for this operator. We then show in Section [3] that this self-adjoint 
differential operator Ls has a discrete spectrum. In Section U] we establish that 
Ls commutes with the operator Ht- This allows us to find the SVD of Ht- In 
Section [5] we then study the accumulation points of the singular values of Ht- In 
particular, we find that Ht is not a compact operator. Finally, we conclude by 
showing numerical examples in Section |6l 
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2. Introducing a differential operator 
Following the procedure in [7] , [8] , we define a differential operator 
Definition 2. 

L{x, dMx) {P{x)i,\x))' + 2(.T - a fi^{x) (2.1) 

where 

4 4 
i=l 1=1 

The four points are all regular singular, and in a neighborhood of each 
the functions (a; — a^) • P'{x)/P{x) and (a: — a^)^ • 2(a; — a)'^/P{x) are complex 
analytic. The term regular singular point is standard in the study of S-L equations, 
see e.g. [16] for this and other terminology and basic properties of S-L equations. 
Consequently, by the method of Fuchs-Frobenius it follows that for A g C any 
solution of L^j — \ip is either bounded or of logarithmic singularity close to any 
of the points a;, see [16] . Away from the singular points the analyticity of the 
solutions follows from the analyticity of the coefficients of the differential operator 
L. More precisely, in a left and a right neighborhood of each regular singular point 
ai , there exist two linearly independent solutions of the form 

oo 

i^i(a;) = |a;-a,ri^5„(a;-a,)", (2.3) 



oo 

i^2{x) = |a;-a,r^^d„(a;-a,)" + fcln|x-a,|i^i(a;), (2.4) 

n=0 



where without loss of generality we can assume — df) = 1. The exponents ai and 
a2 are the solutions of the indicial equation 

+ (Po - l)a + go = 0, 

where 

Po = lim {x - a,)P'{x)/P{x), (2.5) 

qo = lim {x - a,)^2[{x - a f - \]/P{x). (2.6) 

With our choice of P, this gives ai = a2 —0 which implies k ^ 0. For the bounded 
solution in (|2.3p . ai = results in V'i(aj) 0- The radius of convergence of the 
series in (12. 3|) and ()2.4|) is the distance to the closest singular point different from 
ai. In a left and in a right neighborhood of a^, the general form of the solutions of 
(L - A)V' is 

OO 

il,i{x) = loY.^r.{x - a,r (2.7) 

oo oo 

Mx) =hY,dn{x- aiT + £2 In Ix - a,| ^ 6„(x - a,)" (2.8) 

for some constants Ij. Hence we have one degree of freedom for the bounded 
solution, and two - for the unbounded solution. Clearly, for the bounded solutions 
(|2.7p . the coefficients &„ are the same on both sides of a^, since we have assumed 
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60 = 1. However, the bounded part of the unbounded solutions (|2.8p may have 
different coefficients and d+ to the left and to the right of respectively. 

2.1. The Maximal and Minimal Domains and Self-Adjoint Realizations. 

Since we are interested in a differential operator that commutes (on some set to be 
defined) with Ht, we want to consider L on the interval (02, 04). Due to the regular 
singular point 03 in the interior of the interval, standard techniques for singular 
S-L problems are not applicable. It is crucial for our application that we identify 
a commuting self-adjoint operator, for which the spectral theorem can be applied. 
We therefore wish to study all self-adjoint realizations; we follow the treatment 
in Chapter 13 in [21] which gives a characterization of all self-adjoint realizations 
for two-interval problems, of which problems with an interior singular point are a 
special case. 

First of all, one needs to define the maximal and minimal domains on Ij — 
(aj,aj+i) (see Chapter 9 in [21]). Let ACioc{I) be the set of all functions that are 
absolutely continuous on all compact subintervals of the open interval /. Then, 

D„ {^p-.Ij^C: ^, Pi,' e ACiocilj); ^, L^, e £'(/,)}, (2.9) 

-Dj.min := {■'P e -Dj.max : SUpp C (flj, flj + i)}, (2.10) 

and the related maximal and minimal operators are defined as follows: 

^j,max ■ — -^(-^^.max) ■ -^^.max ^ ^ (-^i)' (^-f f ) 

We shall follow essentially the procedure in Chapter 13 in [21], to which we refer for 
more details. On (02, 04), the maximal and minimal domains and the corresponding 
operators are defined as the direct sums: 

Definition 3. The maximal and minimal domains -Dmax; -Dmin C £^(a2,a4) and 
the operators Lmax, ^min are defined as 

-^max • — -^2, max ^" -^3, max? -^min • — -^2,min ^" -^3,min 

-^max • — -^2, max ~t- -^3,max5 -^min ■ — -^2,min ~t" -^3,min 



and therefore 



imax : i'max 'C^(a2, 04), (2-13) 
imin : £'min -> C?{a2, Cli)- (2-14) 



The operator L^mi is a closed, symmetric, densely defined operator in £ (02, 04) 
and Lmax, ^min form an adjoint pair, i.e. L*^^^ = L,nin and L*^;^ = imax- In order 
to define a self-adjoint extension of imin, we need to introduce the notion of the 
Lagrange sesquilinear form: 

[u,v\:^ uPv' -vPu' , (2.15) 

where, at the singular points, 

[u,w](a+) := lim [u,w](a), (2.16) 

[w,w](a~) := lim [u,w](a). (2.17) 
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These limits exist and arc finite for all u, w G Dmax- If we choose u, v d -Dmax such 
that [u,v]{ai) = 1 for all the singular points (a^, , a'^ , Q^r); then the extension 
of -/^min defined by the following conditions 

[^,u]{at)=0 = [^b,u]ia^) (2.18) 

[i^,u]{a^) = [^P,u]{a+) (2.19) 

[V',«](a3) = [V,«](4) (2-20) 
is self-adjoint. We refer to (12. 18^ as boundary conditions, and to (|2.19l) and (|2.20p 
- as transmission conditions. The latter connect the two subintervals (02,03) and 
(03,04). Motivated by the conditions mentioned in Section [1] we define a self- 
adjoint extension of L^in- 

Lemma 1. The extension Ls : D{Ls) — )• £^(02,04) of Ln^m to the domain 
D{Ls) := {ip e Z^max : [iJ,u]{a+) = [tP,u]{a^) ^ 0, 

[^P,u]{a^) = [yj,u]{at),[yj,v]{a^) - [^,«](o+)} (2.21) 
with the following choice of maximal domain functions u,v € D-c^ay. 

u{y) := 1, (2.22) 

«(y):=E n ln|y-a.|, (2-23) 



i=l j^i 

i6{l,---,4} 



is self- adjoint. 



This choice of maximal domain functions gives [u,v\{ai) — 1 for i ~ 1,...,4. 
The boundary conditions simplify to 

P{y)i^\v) = lim P[vW{v) = 0. (2.24) 

For an eigenfunction ip of Ls this is equivalent to ip being bounded at 02 and 
04 (because the only possible singularity is of logarithmic type). Let 0i and 02 
be the restrictions of ip to the intervals (02,03) and (03,04), respectively. Since 
ip is an eigenfunction, on the corresponding intervals 0i and 02 are of the form 
4'i{y) — 4>ii{y) + 4>i2{y) in \y — 03]. Here, the functions (pij are analytic on (02, 04). 
Having this, the transmission conditions can be simplified as follows: 

[V',u](aJ) = [i',u]ia3) 
lim P{y)ip'{y) ~ lim P{y)ij/{y) 

lim 012 (y) = lim 022 (y). (2.25) 

The condition involving v yields 

[i',v]{a-) = [iP,v]{a+) 
lim [iPiy) - v{y){PiP')iy)] = lim [^(y) - v{y){Pil/){y)] (2.26) 

lim 011 (y) = lim 021 (y) (2.27) 

Note that on each side of p.26p the logarithmic terms in 0^2 cancel because of the 
choice of the constants in v. The properties (|2.24p . (|2.25p and (|2.27p are the same 
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as the ones found for /„ in Section [TJ Thus, we have constructed an operator Ls 
for which close to the points 02, as and a4, the eigenfunctions behave in the same 
way that is expected for the /n's. 

Close to 03, an eigenfunction %p is given by 



hi E d^iv-azY 

m=0 



^2iln|y-a3| bmiy-azY 

m=0 



^12 E d+{y 

m=0 



as) 



m=0 



as) 



y <as 
y> as 



d+ 



1 and &o 



(2.28) 

1. The transmission 



where similarly to p.4p . we assume cJq 
conditions require that 

^11=^12, (2.29) 

£21=^2. (2.30) 

We can thus express ip in a neighborhood (03 — 6, 03 + 5), for e > sufficiently small, 
as 



00 

i^{y) = hi + ^21 In \y - as\ ^ b,n{y - as)"" + 

m=0 



et.iy-asY 



(2.31) 



where £^ stands for ^+ = £11^^, when y > as and for £^ ~ ^nd^, when y < as- 

3. The spectrum of Ls 

In order to prove that the spectrum of the differential self-adjoint operator Ls 
introduced in Lemma[T]is discrete, we need to show that for some z in the resolvent 
set, {Ls — zl)~^ is a compact operator. To do so, it is sufficient to prove that the 
Green's function G of Ls — zl, which for z in the resolvent set exists and is unique, 
is a function in £^((02,04)^). This would allow us to conclude that the integral 
operator Tq with G as its integral kernel is a compact operator from £^(02,04) to 
£^(02,04), where Tq is equivalent to the inversion of Ls — zl. 



Lemma 2. The Green's function G{x, ^) associated with Ls — i is in C'^{{a2, (14)^) 
and consequently, {Ls — i)~^ : £^(02,04) — > D{Ls) C £^(02,04) is a compact 
operator. 

Proof. The self-adjointness of Ls is equivalent to Ls — i being one-to-one and 
onto (Theorem VIII. 3 in pTO]). Moreover, the aYs are limit-circle points and thus, 
the deficiency index d equals 4 (Theorem 13.3.1 in [21]). This means that if we 
do not impose boundary or transmission conditions, there are two linearly inde- 
pendent solutions of {L — i)pj = on (02,03) as well as two linearly indepen- 
dent solutions of {L — i)qj — on (as, 04). Note that neither pj nor Qj can be 
bounded at both endpoints because i is not an eigenvalue of the self-adjoint oper- 
ator Lj^s ■■ D{Lj,s) C^{Ij) with D{L^,s) = {V- e I?j.max : lim,„^a+ P{yW{y) = 
lim - P{y)Y''{y) = 0}. By taking appropriate combinations, if necessary, we 

can eliminate the logarithmic singularity at 02 of one of the solutions, and at 04 - 
of another solution. We can thus assume that 

- on (02,03): pi is bounded at 02 and logarithmic at 03, p2 is logarithmic at 
both endpoints; 
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- on (03,04): qi is logarithmic at 03 and bounded at 04, q2 is logarithmic at 
both endpoints. 

We next check the restrictions imposed by the transmission conditions at 03. Close 
to 03, both functions pi and (72 are of the form (I2.8p . Let £11, £21 denote the free 
parameters in the expression for pi and £12, £22 the ones in q2. These can be chosen 
such that they satisfy (j2.29|) and ()2.30|) . Thus, there exists a solution hi{x) on 
(02,04) given by 

Pi{x) for X e (02,03) 



hi{x) 



q2{x) for a; e (03,04) 



that is bounded at 02 and logarithmic at 04. In addition, it is of the form (j2.3ip 
close to 03, i.e. it is logarithmic at 03 and satisfies the transmission conditions 
(|2.25p . (|2.27p there. Similarly, with p2 and qi we can obtain a solution /12 on 
(02,04) that satisfies the transmission conditions at 03 and is of In-ln-bounded- 
type. Thus, imposing only the transmission conditions, we obtain two linearly 
independent solutions of (L — i)h on (02,03) U (03,04). One of them, hi, is 
of a bounded- In-ln- type, and the other one, /i2, is of a In-ln-bounded-type, at the 
points 02, 03, 04, respectively. We are now in a position to consider the Green's 
function G(x,^) of Ls — i- Close to 03, we can write the two functions as hj{x) = 
hji{x) + In |a; — az\hj2{x) with continuous functions hji and hj2. By rescaling if 
necessary, we can assume hi2{a^) = h22{o-3)- We construct G from hi and /12 as 
follows: 

Gix,0 = l ''SJ^Si (3.1) 
\ C2{0h2{x) for x > ^ ^ ' 

where ^ € (02, 03) U (03, 04) and the functions ci(^) and C2(^) are chosen such that 
G is continuous at x = ^ and dG/dx has a jump discontinuity of 1/P(^) at j: = ^: 

ci(0/ii(0-C2(e)/i2(0=0, (3.2) 

ci{0hm-C2{S.)h'2{i)^~^y (3.3) 

In other words, G is the solution of {L — i)G — 5, where S is the Dirac delta function. 
For ^ away from 03, G{x, ^) is continuous in ^ but with logarithmic singularities at 
02 and 04. This can be seen as follows. Consider ^ close to 02. There, we can write 

/i2 (0-^21 (0+^22(0 In I? -02! 

and, since hi is bounded close to 02, it is of the form (|2.7p . i.e. ^.1(02) 7^ 0. Let 
Whi,h2 denote the Wronskian of hi and /i2, i.e. W^/ii,'i2 — ^1^2 ~ h[h2- For ci and 
C2 we obtain 

The denominator in the above expressions is bounded by 

p{mim'2{o - h2m'i{0) = oa^ - 02) in - 02I) + /ii(o^22(Op(o, 
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where p{^) — P{£,)/{£, ~ 02) and hi{a2)h22{a2)p{a2) 7^ 0. Thus, m a neighborhood 

of 02, 

ci(0=O(ln 1^-02!), (3.6) 

C2(0 = 0{1). (3.7) 

Similarly, since /i2(a4) 7^ 0, close to 

ci(0=O(l), (3.8) 

C2(e) -0(ln|?-a4|). (3.9) 

For each fixed ^ G (a2,a3) U (03,04), G{x,^) as a function in x is continuous on 
[02, 03) U (03, 04] and has a logarithmic singularity at 03, due to the singularities in 
hi{x) and h2{x). It remains to check what happens as ^ — 03. We need to make 
sure that the functions ci(^) and C2(^) behave in such a way that G G £^((02, 04)^). 
Therefore, we derive the asymptotics of ci(^) and C2(^) as ^ Og". For ^ = 03 — e 
and small e > 0, equation (13. 2p becomes 

Ci(a3 - e)hi{a3 - e) - 02(03 - e)ft.2(a3 - e) = 0. 

Since close to 03, hi — hn + h-a ln(e) and the hij are continuous, the ratio ci/c2 is 
of the form 

a + 61n(e) 
c + dln(e)' 

where h and d are non-zero (because the logarithmic singularity is present). Thus, 
the ratio tends to the finite limit b/d as, e ^ Q. Conditions (|3.2p and p.Sp together 
imply: 

ft,i(o3 - e) 



C2(a3 - e) 



P(o3 - e)M^/ii,/i2(a3 - e) 
/til (as - e) + ^12 (03 - e) ln(e) 
ri(e) + e • r2(e) 



where 



r'i(e) = -ip(o3 - e){h2ihi2 - /i22^ii)(a3 - e), 
'^2(e) = 0(l), 

and ^.12(03) 7^ 0. If ri(0) 7^ 0, then C2 is of order C'(ln(e)), removing a possible 
obstruction to square integrability of G. 
Suppose ri(0) = 0, i.e. 

/i2i(a3)^i2(a3) - ^22(03)^11(03) = 0. 

This would imply 

/iii(a3) -C-/i2i(a3), (3.10) 

/ii2(a3) -C-/(,22(a3) (3.11) 

for some constant C. By assumption, /ii2(o3) = h22{aj,), so that C = 1. Now if 
both (|3?T0l) and ([3lT|l hold for C = 1, the function defined by 

hi{x) for X e (02, 03) 
h2{x) for X e (03, 04) 

would be a non-trivial solution of (L5 — i)h — (fulfilling both boundary and 
transmission conditions), i.e. i would be an eigenvalue of L5. But this contradicts 



h{x) 
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the self-adjointncss of Ls- We can thus conclude that ri(0) ^ 0. This shows that 
02(03 — e) is of order C'(ln(e)) and therefore also C2 • = ci = C'(ln(e)). 
Analogously, we can find the same asymptotics of ci(^) and C2('C) as ^ — )■ a^. 
Therefore, the properties of the Green's function G{x,£^) can be summarized as 
follows: 

- G(-,^) has logarithmic singularities at a2, and 04, 

- G{x,-) is of logarithmic singularity at as, 

- away from these singularities G'(x,^) is continuous in x and ^. 

Thus, G is in £^((02,04)^). Hence, Tq ■ £^(02,04) — )■ £^(02,04) is a compact 
Fredholm integral operator. □ 

From this we conclude: 

Proposition 1. The operator Ls has only a discrete spectrum, and the associated 
eigenfunctions are complete in £^(02,04). 

Proof. By Theorem VIII. 3 in [15], the self-adjointness of Ls implies that for the 
operator [Ls — i) : D{Ls) — >■ £^(02,04) we have 

Ker(Ls - = {0}, (3.12) 

Ran(_L5 — i) = £■^(02, 04). (3.13) 

Consequently, {Ls ~ i)^^ : £^(02,04) — > D{Ls) is one-to-one and onto. Moreover, 
it is a normal compact operator and thus we get the spectral representation 

00 

{Ls-l)-'f = Y,\n{fJn)fn, (3.14) 

where {/n}neN is a complete orthonormal system in £^(02,04). This can be trans- 
formed into the spectral representation for L5: 

is/ = E(ir + (3.15) 

□ 

Clearly, the eigenfunctions /„ of Ls can be chosen to be real-valued. The com- 
pleteness of {/„}„gN is essential for finding the SVD of Ht- Another property that 
will be needed for the SVD is that the spectrum of Ls is simple, i.e. that each 
eigenvalue has multiplicity 1. 

Proposition 2. The spectrum of Ls is simple. 

Proof. From the compactness of (L5 — i)~^ , we know that each eigenvalue has 
finite multiplicity. Suppose /i and /2 are linearly independent eigenfunctions of Ls 
corresponding to the same eigenvalue A e M. Then, on all of (02,03) U (03,04) the 
following holds 

fiLf2 - fiLfi = 0. (3.16) 

Consequently, 

= fM - f2Lfi = fiiPf^y - f2{pf[)' 

= [/l,/2]'. 
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Thus, [/i,/2] is constant on both (02,03) and (03,04). From the boundary con- 
ditions that /i and /2 satisfy, we find that [/i,/2](a2) = = [/i,/2](o4), which 
imphes [/i,/2] = on (02,03) U (03,04). Since [/i,/2] = P{f[f2 - /1/2), we get 
that 

f[f2 - /1/2 = on (02, 03) U (03, 04). (3.17) 

The functions /i and /2 satisfy the transmission conditions at 03. Consequently, 
they can be written as 

fi{x) = fii{x) +/i2(a;)ln|a;- 03I, 

f2{x) = f2iix) + /22(a;)ln |a; - 03], 

in a neighborhood of 03, where are continuous. Since the one-sided derivatives 
f-j are bounded at 03, equation (I3.17P imphes 

(/l2/21-/ll/22)(x) ^ 1^ _ ^^1^ ^ ^ 

X - 03 

Note that the terms containing hi \x — a3|/(a; — 03) canceL Taking the hmit a; — >■ 03 
in p.lSp . we obtain 

/l2(03)/2l(03) - /ll (03)722(03) = 0. 

Thus, for some constant C: 

/ii(a3) \ ^ Q f /2i(a3) 



712(03) / \ 722(03) 

If we take /i on (02,03), then /ii(o3) and 712(03) define a singular initial value 
problem on (03,04) that is uniquely solvable (Theorem 8.4.1 in [H]). Thus, 7i = 
C ■ 72 on (03, 04). Now, on the other hand, by considering 7i on (03, 04), the values 
711(03) and 712(03) define a singular initial value problem on (02,03) which has a 
unique solution. Hence, fi = C ■ f2 on (02,03) U (03,04) in contradiction to our 
assumption. □ 

4. Singular value decomposition of Ht 

Having introduced the differential operator Ls, we now want to relate it to 
the truncated Hilbert transform Ht- The main result of this section is that the 
eigenfunctions of L5 fully determine the two families of singular functions of Ht- 
We start by stating the following 

Proposition 3. On the set of eigenfunctions {fn}neN of Ls, the following com- 
mutation relation holds: 

{HTL{y,dy)fn)ix) = L{x,dx){HTfn){x) 7or a; e (oi, 02) U (02, 03). (4.1) 

Sketch of proof. This proof follows the same general idea as the proof of Proposition 
2.1 in 0. We therefore provide full details only for those steps where additional care 
needs to be taken because of the singularity at 03. The steps that are completely 
analogous to those in the proof of Proposition 2.1 in are only sketched here. 

Let ■0 € {7n}nGN- The boundedness of ^ at 02 and 04 implies that PV' ~^ ^-ud 
Pip there. Moreover, the transmission conditions at 03 guarantee that Pip' is 
continuous at 03. With these properties, the commutation relation for x G (01,02), 
i.e. where the Hilbert kernel is not singular, can be shown similarly to the proof of 
Proposition 2.1 in [7]. 
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Next, let X € (02,03). The main difference from the proof of Proposition 2.1 
in [8] is that now the eigenfunctions are not in C°°([a2, 04]), but are singular at 
03. However, the fact that we exclude the point x = ag, allows us to always have a 
neighborhood of x away from 03 on which ip is bounded. We further note that ip € 
C°°{[a2, 03) U (03, 04]). Since the Hilbert kernel is singular, we need to use principal 
value integration and introduce the following notation: Ie{x) := [02, x— e]U[a;+e, 04]. 
Here e > is so small that {x — e,x + 2e) C (02, 03), i.e. the e-neighborhood of x 
is well separated from 03. Then, 

-{PiyWiy))' , 2(y-^)2^(y)i 



■K{HTL{y,dy)ip){x) 



lim 

e->-0+ 



y 



dy. 



For the first term under the integral, we integrate by parts twice and plug in the 
boundary conditions. Again, we use that Pip' — )■ and Pip — >■ at 02 and 04; 



dy 



{Pij'){x -e) + {PtP'){x + e) {Ptp){x - e) - {Pi^){x + e) 



h[x) 



^ 2P{y)-P'{y){y-x) ^ 

VXy) 7 X3 dy. 

(y - xy 



(4.2) 



The integral on the right-hand side of (|4.2p can be related to the derivatives of 
/ i'{y)/iy ~ x)dy. In |5] similar relations (cf. eq. (2.7)) were obtained from the 
Leibniz integral rule, using explicitly that the integrand was continuous. In our 
case, the function ip is no longer continuous because of the singularity at 03. We 
can generalize the argument of [5] by invoking the dominated convergence theorem 
and rewrite the last term in (14.21) as follows: 



^^yfP(y)-ny)iy-x),^^ 



I,{x) 



P{x) 
+ P'{x) 



{y - xY 

^ I ^dy 



tp'{x — e) + ^'{x + e) 1/1(2; — e) — ^j{x + e) 



dy 



^{x — e) + ip{x + e) 



dx'^ J y — X 

I,{x) 

d_ r ^Pjy) 
dx J y- X 

h(x) 

„ , , Ay - aY - ix - (tY , 

my)- — - — ^ -dy 

y-x 

h{x) 

Putting all pieces together, we obtain: 

TT{HTL{y,dy)ip){x) = 

{P^')ix - e) + {Pi;'){x + e) {P^)ix - e) - {P^){x + e) 



lim \ 

e-)-0+ L 



e 

ip'{x — e) + ip'{x + e) ijj{x — e) — ipix + e)n 



P{x)[ 

c 

P'(x)^^fc^l±^^(^+L(.,d.) 



h{x) 



y-x 



dy} 



SPECTRAL ANALYSIS OF THE TRUNCATED HILBERT TRANSFORM 



15 



The eigenfunction tp is in C°°[a2,x + 2e). Following [8], we can thus express the 
boundary terms in the above equation by Taylor expansions around x and make 
use of the fact that the boundary terms consist only of odd functions in e. The 
boundary terms are then of the order 0{e). We thus have 

TT{HTL{y,dy)ij){x) = \im L{x,d,) f ^^dy. (4.3) 

<^^o+ Jh{x) y-x 

Since for e > sufficiently small, ■0 G C°°{\x — e,a; + e]), one can interchange the 
limit with L{x,dx) as in [8]. □ 

Because the spectrum of L5 is purely discrete, we have thus found an orthonor- 
mal basis (the eigenfunctions of Lg) {/„}„gN of £^(02,04) for which (|4.ip holds. 
Let us define g„ :— HTfn/\\HTfn\\c^(ai,a3)- Then, in order to obtain the SVD 
for Ht (with singular functions /„ and (?„), it is sufficient to prove that the g^s 
form an orthonormal system of £^(01,03) (they will then consequently form an 
orthonormal basis of £^{ai,a3), see Proposition [5]) . 

The orthogonality of the g^s will follow from the commutation relation. Since 
/„ is an eigenfunction of L5 for some eigenvalue A„, we obtain 

L{x,dx)gn{x) = A„.g„(a;), x G (01,02) U (02,03). 

Similarly to L5, we define a new self-adjoint operator that acts on functions sup- 
ported on [oi, 03]: 

Definition 4. Let Z^max := £'i,max + £)2,inax and Lmin := ii,min + i2,min- The 
operator Ls : D{Ls) £^(01,03) is defined as the self-adjoint extension of L^-m, 
where 

D{Ls) := {V' e £'max : [^,u]{a+) = [7/',m](o~) = 0, 

[i>,u]{a^) = [i,,u]{a+),['^,v]{a^) = [?/',f](o+)} (4.4) 

with the maximal domain functions u,v G -Dmax os in ()2.22p . (I2.23p . 

The intuition then is the following. The function /„ is bounded at 02 and 
logarithmic at 03, where it satisfies the transmission conditions. Consequently, 
as will be shown below, gn is bounded at 03, logarithmic at 02 and satisfies the 
corresponding transmission conditions at 02. Clearly, it is also bounded at oi. 
Thus, gn is an eigenfunction of the self-adjoint operator L5. As a consequence, the 
(7„'s form an orthonormal system. 

Proposition 4. If Lsfn = Kfn, then gn ■= Hrfn/WHrfuWc'^ (01,03) an eigen- 
function of Ls corresponding to the same eigenvalue 

Lsg-n = A„g„. (4.5) 

Proof. First of all, the commutation relation for /„ yields 

L{x,d^){HTfn){x) - {HTL{y,dy)fn){x), 

L{x,dx)gnix) = Xngn{x), X G (01,02) U (02,03). 

What remains to be shown is that gn satisfies the boundary and transmission 
conditions. Therefore, we consider p.v. J^* fn{y)/ {y — x)dy for x close to oi, 02 
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and 03. In a neighborhood of ai away from [a2, 04] this function is clearly analytic. 
Now, let X G (02 — e, 02 + e) and choose ft, such that 02 + e < 03 — ft. Then, 

p,,fMyldy = p.v. [ ^dy+ f i^dy. (4.6) 

J y~ X J y- X J y-x 

The second expression on the right-hand side is analytic on x G (a2 — e, 02 + e). For 
the first term we can write 

as — h a^—h as — h 

f fn{y) , f fn{y)-fn{x) . . f 1 , 

p.v. / dy^p.v. / dy + }n{x) ■ p.v. / dy 

J y- X J y-x J y- X 

0.3— h 

r f^(y) — fn{x) 

^p.v. / — dy + /„(a::)(ln |a3 - ft - xl - ln|a2 - a;|). 

J y-x 

0-2 

The first expression on the right-hand side is bounded, since /„ is Lipschitz con- 
tinuous on [02,03 — ft.]. Moreover, it is also continuous, which follows from the 
smoothness of {fn{y) ~ fn{x))/{y — x). Thus, close to 02, gn has a logarithmic 
singularity of type 

9n{x) = gn,l{x) + gn.2{x) In \x - a2\, 

where g„ i are continuous at 02 , thus gn satisfies the transmission conditions (j2.25l) , 

It remains to check the behavior of gn close to 03. Let x = — S for a small 
5 > 0. As before, for a fixed ft > 0, we can split the integral 



p.v. 



a4 as'—h a4 as+Zi 

l^dy^p.v.f I + I )Myldy+p.v. I l^dy 

y-x \ J J J y-x J y-x 

a2 <i2 a3+/i as—h 



where the first term on the right-hand side is analytic for S < h. Thus, we are left 
with the second expression on the right-hand side. We can choose ft > 0, such that 
/„ is of the form (j2.3ip in (03 — ft, 03 + ft). In this neighborhood around 03, we can 
express /„ as 

fn{y) = fn,l{y) +/n,2(2/)ln|?;~ 03], 

where both /„_i and fn,2 are Lipschitz continuous. Then, the integrals 

fn.iiy) 



p.v. / '-^^^dy (4.7) 

J y-x 

a3+h 

p.v. I fnAy)^^\y--^\^y (48) 
J y-x 

a3~h 

are bounded as x tends to 03. This follows from the results in [5], Sections 8.1 and 
8.5. More precisely, by taking the integrals (|4.7p and (|4.8|) over [03 — ft, 03] and 
[03, 03 -|- ft] separately, one obtains that the singularities at 03 cancel. □ 

Since the spectrum of is simple, we can conclude that the gnS form an 
orthonormal system and thus the following holds: 
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Theorem 1. The eAgenf unctions /„ of Ls, together with gn '■— HTfn/\\HTfn\\c^{ai,a3) 
and an ■= \\HTfn\\c^{ai,a3) form the singular value decomposition for Ht'- 

Hrfn = (Jngn, (4.9) 
H^gn = CTri/n- (4-10) 

5. Accumulation points of the singular values of Ht 

The main result of this section is that and 1 are accumulation points of the 
singular values of Hx- To find this, we first analyze the nuUspace and range of Ht, 
which will also prove the ill-posedness of the inversion of Ht- First, we need to 
state the following 

Lemma 3. // the Hilhert transform of a compactly supported f G C^{a, b) vanishes 
on an open interval {c,d) disjoint from the object support, then f — on all o/M. 

Sketch of proof . A similar statement (and proof) can be found in 1 . The main 
difference is that here we consider a more general class of functions /. By dominated 
convergence, / £ L^{a,b) implies that for any z E ft = C\((— oo,a) U (fo, oo)), the 
function g(z) — f{x)/{x — z)dx is differentiable in a neighborhood of z. Thus, 
g is analytic on i7. The statement then follows in the same way as Lemma 2.1 in 

m- □ 

With this property of the Hilbert transform, we can obtain results on the nuUspace 
and the range of Ht'- 

Proposition 5. The operator Ht '- £^(02,04) — > £^(01,03) has a trivial nuUspace 
and dense range that is not all 0/ £^(01,03), i.e. 

Ker{HT) = {0}, (5.1) 

Ran{HT)^ £^{01,03), (5.2) 



Ran{HT) =£^(01,03). (5.3) 
Proof of dsn). Suppose iJr/ = 0. Then 

Hxia2Mi]f = on (01,02), 

and by Lemma [31 / = on all of [02,04]. Thus, / g £^(02,04) can always be 
uniquely determined from HTf- 

Proof of ()5.2|) . Take any g € £^(01,03) that vanishes on (01,02) and such that 
Il5li£2(aj_a3) ^ 0. Suppose g G Ran(iJT). By Lemma [31 if / G £^(02,04) and 
HtJ = g, then / is zero on [02,04]. This implies that 5 = on (01,03), which 
contradicts the assumption ||g|| 7^ 0. 

Proof of (|5.3p . The operator H^ is also a truncated Hilbert transform with the 
same general properties. By the above argument, Ker{H^) = {0}. Thus, 
Ran{HT)-^ = {0}. □ 

Equation (|5.2p shows the ill-posedness of the problem. It is not true that for every 
g G £^(01,03) there is a solution / to the equation i?T/ — 9- Since Ran(7JT) is 
dense, the solution need not depend continuously on the data. Thus, our problem 
violates two properties of Hadamard's well-posedness criteria [3]. These are the 
existence of solutions for all data and the continuous dependence of the solution on 
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the data. We now turn to the spectrum of H^Ht- In what fonows, || • || denotes 
the norm associated with £^(M), and (•, •) denotes the £^(M) inner product. We 
begin with proving the fohowing 

Lemma 4. The operator H^Ht has norm equal to 1. 
Proof. From \\H\\ = 1, we know that \\H^Ht\\ < 1. Since 

\\H^Ht\\ ^ sup \\H^HtM\, 
IIV'll=i 

finding a sequence ■(/'„ with WipnW — 1 and \\H^HT^n\\ 1 would prove the asser- 
tion. 

Take a compactly supported fimction G £^([—1,1]) with = 1 and two 
vanishing moments, ip{x)dx — = J ^^ x-^j{x)dx. From this, we define a family 
of functions, such that the norm is preserved but the supports decrease. More 
precisely, for a > 2/(03 — 0.2), we set 

^,(x) = V^^(a(x-^^i^)). (5.4) 

These functions satisfy W^/jaW = 1 and supp Va = [^^4^ - ^, + ^1 ^ 03]- 

For their Hilbert transforms we obtain 

We can write 

HrHxi^a = -Xla2,ai]HX[ai,a3]HXla2,a4]i'a 

= ~X[a2M4]H{I -{I- XlaiM3]))HXla2.a4]'4^a 
= V-a + X[a2,ai]H(I - X[ai .ag] )^^X[a2 ,04] "^a ; 
(/ - H^HT)iJa = -X[a2,a^]H{I ~ X[ai ,03] )i^X[a2 ,04] V'a ■ (5.6) 

Consider the £-^-norm of the last expression 

II (/ - H^HT)iJa\\'' = ||X[a.,a4]^^a " X[a, ,a3] )^^X[a. ,a4] V'a f 
< - Xlar,a,])Hxla2,a^]'^aV 

\{Hi;a){x)\^dx 

( — oo,ai)U(a3,oo) 

\{Hi,){a{x-^^)\^dx 

( — oo,ai)u(a3,oo) ^ 

\H^P\^dy+ / \Hi,\^dy. (5.7) 

-oo Ja-(a3- °^2°^ ) 

Because of the ordering of the a^'s, we have that ai — < and 03 — 22±£3 > q. 

Since ip has two vanishing moments, H^p asymptotically behaves like l/|j/P and 
hence, both integrals in (|5.7p are of the order 0{a~^). Thus, given any e > 0, one 
can find a > 2/(03 ~ 02) such that 

Wii - h^Ht)M < ^- 

Consequently, 
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Therefore 

WH^HTi^aW 1 as a -H- oo, 
which imphes that \\H^Ht\\ =1- □ 

We are now in a position to prove 

Theorem 2. The values and 1 are accumulation points of the singular values of 

Proof. First of all, and 1 are both elements of the spectrum a{H^HT)- For the 
value 0, this follows from Kan{H^HT) C Ran(iJ^) ^ £^(02,04). Moreover, since 
||ffyi?T|| = 1 and H^Ht is self-adjoint, the spectral radius is equal to 1. Thus, 
1 G a{H*HT). 

The second step is to show that and 1 arc not eigenvalues of H^Ht- 

is not an eigenvalue: If H^Hrf = 0, then WHrfW^ = {f,H^HTf) — 0. Since 
KeriHr) = 0, this implies / = 0. Thus, KeiiH^Hr) = {0}. 

1 is not an eigenvalue: Suppose there exists a non- vanishing function / e £^(02, 04), 
such that 

X[a2 ,04] / = /• 

Then, 

\\HX[a2,ai]ff = Wff -{X[a2.ai]HX[aua3]HX[a2.ai]f: f) 
(X[ai ,03] ^X[a2 ,a4] -^X[a2.a4]f} 
IIX[ai,a3]-^^X[a2,a4]/|| ■ 

This implies that Hx[a2,ai]f is identically zero outside [ai,a3]. By Lemma[3l this 
implies / = X[a2,ai]f — 0, contradicting the assumption / ^ 0. Therefore, and 
1 are accumulation points of the eigenvalues of H^Ht and consequently, of the 
singular values of Ht- □ 

Since the singular values of Ht also accumulate at a point other than zero, the 
operator Ht is not compact. 

6. Numerical illustration 

We want to illustrate the properties of the truncated Hilbert transform with 
overlap obtained above for a specific configuration. We choose oi = 0, 02 = 1.5, 
as = 6 and 04 = 7.5. First, we consider two different discretizations of Ht and 
calculate the corresponding singular values. We choose the first discretization to be 
a uniform sampling with 601 partition points in each of the two intervals [0, 6] and 
[1.5,7.5]. Let vectors X and Y denote the partition points of [0,6] and [1.5,7.5] 
respectively. To overcome the singularity of the Hilbert kernel the vector X is 
shifted by half of the sample size. The z-th components of the two vectors X and 
Y are given by Xi = -^{i + \) and Yi = 1.5 + y^z; Ht is then discretized as 



{HT)i,j = (l/7r)(Xi — Yj), i,j — 0,... ,600. Figure 5(a) shows the singular values 
for the uniform discretization. We see a very sharp transition from 1 to 0. 

The second discretization uses orthonormal wavelets with two vanishing mo- 
ments. Let (j) denote the scaling function. For the discretization we define a finest 
scale J — —7. The scaling functions on [1.5,7.5] are taken to be 4>-7,k for integers 
k = 192,..., 957, i.e. such that supp 4'-7,k C [1.5,7.5]. On the interval [0,6] the 
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scaling functions are shifted in the sense that we take them to be for inte- 

gers £ = 0, . . . , 765, i.e. such that supp (/)_y C [0, 6]. Figure [5(b)] shows a plot 
of the singular values of this wavelet discretization of Ht- Although the transition 
is not as sharp as in 5(a) the singular values in both cases very clearly accumulate 
at and 1. 



300 400 500 fiCK) 



100 200 3m 400 500 



(a) Uniform discretization of size 601 X 601. (b) Wavelet discretization of size 766 X 766. 



Figure 5. ai = 0, 02 = 1.5, 03 = 6, 04 = 7.5. Singular values of 
two discretizations of Ht ■ 

Next, we consider the singular functions. Figure [5] shows the singular functions 
of the uniform discretization for singular values in the transmission region between 
and 1 . Figure [7] illustrates the behavior of singular functions for small singular 
values. As anticipated, they are bounded at the two endpoints and singular at the 
point of truncation. Figure |8] gives two examples of the close to linear behavior in 
a log-linear plot of the singular functions. In agreement with the theory in Section 
m these plots confirm that the singularities are of logarithmic kind. 

Based on the numerical experiments conducted, we make the following obser- 
vations on the behavior of the singular functions and singular values. First, the 
singular functions in Figures [6] and [7] have the property that two functions with 
consecutive indices have the number of zeros differing by 1. Moreover, the zeros are 
located only inside one subinterval Ij. Furthermore, the plots show that singular 
functions with zeros within the overlap region correspond to significant singular 
values, whereas those which have zeros outside the overlap region correspond to 
small singular values. Finally, we remark that singular functions for small singular 
values are concentrated outside the ROI I2 = [02,03]. 
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(a) Singular functions /44g and fag. (b) Singular functions /450 and /451. 




(c) Singular functions 3448 and 9449. (d) Singular functions 3450 and 9451. 



Figure 6. Consecutive singular functions for the uniform dis- 
cretization with 3, 2, 1 and no zeros within the overlap region. 
The corresponding singular values are (7448 — 0.999963, (7449 — 
0.99 8 7 82, (7450 = 0.966 1 92, 0-451 = 0.542071. 
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(a) Singular functions 7452 and 7453 . (b) Singular functions 7454 and 7455 . 




(c) Singular functions (7452 and 9453. (d) Singular functions (7454 and 9455. 



Figure 7. Consecutive singular functions for the uniform dis- 
cretization with 1, 2, 3 and 4 zeros outside the overlap region. 
The corresponding singular values are 0-452 — 6.29189 • 10"'^, 
cr453 = 2.83 5 33 •10"^ a454 = 1.18274- 10-^ (7455 = 4.83357- IQ-i". 

0,00 h 



(a) 



11.050 0.1011 



(b) 



Figure 8. Log-linear plot, i.e. with a logarithmic x scale of 
the singular functions 3459 (left) and 5453 (right) on the interval 
[1.5,2.3]. 
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